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ABSTRACT 

The physical conditions in a high density hydrogen gas heated 
and ionized by supratherraal protons are investigated , with appli- 
cation to the gas in the nuclear region of Seyfert galaxies. The 
gas is assumed optically thick to Lyman and Balraer line radiation. 
J^echanisms by which the radiation from the gas can balance the 
heating by the fast protons are investigated/ and minimum values 
for the mass of gas are estimated. Under certain conditions, the 
suprathermal atoms themselves can cool the ambient gas by rescat- 
tering the line radiation into the optically thin region in the 
v;inqs of the line. This mechanism, which we call optical rever- 
beration, can enhance the broad component of the hydrogen lines 
produced by inelastic atomic collisions and yield line widths con- 
sistent with those observed in Seyfert and quasar spectra. If 
this mechanism is important, the Lyman decrement can depend strongly 
on the temperature and density of the ambient hydrogen. We also 
discuss the possibility of achieving dynamic equilibrium of the 
ambient gas by balancing the momentum transfer from the suprathermals 
with gravitational attraction of a massive central source. 
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I. introduction 

A radially' streaming proton model in which suprathermal pro- 
tons interact with a stationary ambient gas of hydrogen atoms to 
give rise to emission lines, has been successful in understanding 
the broad wings of Balmer lines in Seyfert galaxies (see Ptak and 
Stoner 1973) . In the model, the physical properties of the ambient 
hydrogen are free parameters that have insignificant effect on the 
Lyman and Balmer emission line profiles. Variations on the model 
have produced agreement with the observed emission line profiles 
of other ions in quasars (Stoner, Ptak and Ellis 1974) and demon- 
strated that considerably different geometries can produce nearly 
identical profiles (Stoner and Ptak 1974) . Thus it appears that 
i the essential feature of the radially streaming proton model in 
producing agreement with observed profiles is the role of the in- 
elastic atomic collision processes. 

A flux of fast protons large enough to produce emission wings 

\ 

of the observed intensity by the atomic collision mechanism should 
play a major role in determining the physical state of the ambient 
gas. The kinetic energy and the momentum of the fast protons will 
be transferred to the gas, producing heating, ionization and, 
possibly, large-scale dynamic effects. Here we discuss a simplified 
model that is designed both to assess the effect of the suprathermal 
protons on the temperature and state of ionization of the gas and 
to investigate the ways and foms in which the kinetic energy of 
the protons can eventually appear in the radiated electromagnetic 
spectrum of the object. By using the observed spectra of Seyfert 
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galaxy nuclei and by requiring physical self consistency, we 
estimate the physical properties of the nuclear gas from which 
the broad emission lines originate. 

IX. Description of the Model 

The model we consider assumes uniform physical properties 

for the ambient gas. Detailed self consistent calculations by 

Kimmer (1975) show that, except for surface effects, this is a 

reasonable approximation, even when the velocity distribution of 

the fast protons varies considerably from place to place. The 

suprathermal protons are assumed to be homogeneously distributed 

and moving in random directions; they are characterized by an 

injection energy and an overall rate of injection. While fits 

to observed profiles require some coherence in the directions of 
suprathermal ion motions , the observations give few other clues 

about geometry, and this coherence will affect momentum transfer 

rates rather than heating rates. Thus, geometrical considerations, 

necessary for fitting observed profiles, are disregarded, except 

that the linear dimension and the total volume of the gas are 

parameters. The gas is assumed to be purely atomic hydrogen, 

partially ionized. 

If the ambient gas were optically thin, the emission from 
the ambient atoms would contribute strong central cores to the 
hydrogen lines, with an intensity greater than that contained in 
the broad component. Since this is not observed, we assume that 
the density of the ambient gas, its linear dimension, and the 
n = 2 population are all sufficiently large to make the gas optically 



5 


thick in both Lyma,n ^nd Baliter lines. We J^urther assume that the 
gas is optically thick In the Lyman continuum, hut transparent to 
the Balmer continuum, so that only a part of the energy delivered 
when an atom is ionized can directly escape via recombination. 

Then the line radiation is trapped until it can be degraded into 
an optically thin region by one of the mechanisms we describe below. 

Most of the energy of the suprathermal protons is first 

delivered to the electrons of the gas, with the rest going into 

excitation of the ambient atoms. The electron temperature rises 

until the electrons are able to transfer energy to the ambient gas 

atoms at the same rate as they receive it from the suprathermals . 

Both thermal bremstrahlung and direct heating of dust by the 

electron gas are generally negligible. Radiative cooling through 

forbidden line emission from heavy elements is suppressed by col- 

7 12 -3 

lisional de-excitation at the densities we consider (10 to 10 cm ) 

Next the energy goes by way of the atoms into the trapped 
radiation field in the vicinity of the emission line centers. Since 
each Lyman and Balmer photon will be absorbed and reemitted a very 
large number of times before being degraded or escaping, the atoms 
are very efficient "heat pipes" connecting the heat baths in the 
radiation near the line centers. Thus, we assume the radiation in 
the optically thick region near the line centers rises to the level 
specified by the Planck formula with a radiation temperature that 
also describes the distribution over atomic levels of the ambient 
hydrogen. The free electrons and the Boltzmann distribution over 
atomic levels are described by different temperatures whose relative 



values are to be fixed by requiring that the net rate at which the 
electrons deliver energy to the radiation field via the atoms is 
balanced by the rate at which the radiation field is cooled. 

We consider three ways by which the energy can finally leave 
the gas. One possibility is that the two photon decay of the 
hydrogen 2s state acts as the dominant cooling mechanism, contri- 
buting an ultraviolet continuum to the radiated spectrum. Second, 
if enough dust is present, the line radiation can be absorbed and 
reemitted in the infrared. Finally, the energy may, under certain 
conditions, be absorbed and reemitted by the suprathermal atoms 
themselves : in this case the energy will come out in the broad 
component of the hydrogen emission lines. The proportion of the 
energy radiated in each of these three forms will depend on the 
temperature, the density of ambient gas, the flux of suprathermal 
particles, and the amount of dust. If the third mechanism is the 
dominant mechanism in determining the energy balance, then it can 
also have a large effect on the Lyman and Balmer decrements, and 
the relative hydrogen line intensities will be determined by the 
temperature and density of the ambient gas. 

III. Ionization and Energy Balance 

For a given temperature of the gas of free electrons, the 
ionization fraction is fixed by a balance between ionization and 
recombination rates for the ambient hydrogen: 
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where N, is the ambient number density of La photons , Njj is the 

number density of hydrogen atoms in the n = 2 state , 

the cross section for ionization of n = 2 hydrogen by La photons, 

N„ is the number density of hydrogen (protons plus atoms) , f ^ is 
the ionization fraction of the ambient gas, n^2 

recombination coefficient for hydrogen, and Y (Tg) is the corres- 
ponding rate coefficient for ionization by the free electrons. The 
effect of the suprathermal protons on the ionization rate is repre- 
sented in equation Cll by a term containing a number density of 
suprathermal particles N^, which is related to the overall in- 
jection rate J (i.e., suprathermal protons injected per unit time), 
the time required for an injected proton to thermalize, t, and the 

total volume of ambient gas by 

N = Jt/V (2) 

s 

The slowing dovTn time will depend on the injection energy, 
the density of ambient gas, and the ionization fraction. If 
Coulomb collisions with free electrons is the dominant mechanism 
for decelerating the protons, then the slowing down time will be 


approximately 


3.2 X 1q8 / 

\^25kevJ 


(sec) 


(3) 


Equation (3) is an approximation based on the formula by Butler 
and Buckingham (1962) for the deceleration of a suprathermal charged 
particle in a plasma; it neglects some dependence of the average 
charge of the suprathermal proton on its velocity and the slight 
dependence of the DeBye screening factor on density , temperature , 
and velocity. 
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For th.e remaining- unspecified quantities in equation Cl), we 
use tile following approximations: 

C4) 

(5) 

(6) 

each in units sec ^ cm^. Equation (5) is the ionization rate co- 
efficient of Kalkofen and Strom (1966) and equation (6) represents 
an approximate average over ionization rates due to suprathermal 
atoms and protons. 

We assume the ^lominant mechanism by which the free electron 
gas exchanges energy with the radiation field is the collisional 
©xcitation and de— excitation of the n = 2 state of the hydrogen 
atom. If there were nearly complete thermal equilibrium betv/een 
the electron gas and the radiation field, the excitation and de- 
excitation rates would be in nearly perfect balance, and the two 
temperatures, T and T , would be equal; however, the two-photon 
decay of the hydrogen 2s state provides a minimum rate of cooling 
that reduces the effective temperature of the radiation relative 
to that of the electron gas. If we equate the rate of collisional 
excitation of the n = 2 state with the sum of the rates for colli- 
sional de-excitation and two-photon production, we obtain an appro- 
ximate relationship between the two temperatures: 



c^CTg) = 2.07 X 10"^^C6,1 - .961og3^Q(Tg))/Tg^ 
rCT^i - 5.3 X 10”^^Tj^exp(.-1.58 x lO^T ) 

0 G ^ 

—8 

<cr. Csiv> = 2.0 X 10 ; ca._„(La) = 1.8 x lo 

ion ion 
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where, as a aimplWyi-nS aaswoptian, the cross section for exortatron 

O* n = 2 By electron impact has Been taken to vary Inversely with 

. • 4^ T% 4- "T ^inpt*civ is 3.bov^ tlur 6SllO 1(3- • 

the impact energy -when the electron 

Xf other cooling mechanisms are significant, then the depop- 
nlation of the n = 2 state By emission of Lyman alpha photons will 
not be entirely balanced by the reverse absorption process, and 
the radiation temperature will be correspondingly smaller than that 
given by equation (7) . The remaining requirement for a steady 
state of ionisation and temperature is that the total rate of heating 
of the gas by the suprathermals equals the rate at which radiated 
energy escapes the gas. The overall rate at which energy is de- 
livered to a unit volume of the gas is clearly J E,/v, or, using 
equations (2) and (3), the heating rate is: 


S O 


«S^iVo 
3.2 X 10* 



( 8 ) 


If we assume that the Lya photons produced by excitation of 
the ambient gas are inraediately degraded (through absorption by 
dust, for example), then the above heating rate must be in balance 
with the rate at which energy is transformed into Lya photons by 
electron excitation, which is 




(9) 


Where c,,, (e) is the cross section for oollisional excitation of 
Lva bv electrons, hv„ is the energy of a Lya photon, and the average 
is over the Boltsmann distribution of the electron velocities. 
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Equating expressions C8) and (9) gives us the energy balance, and 
this, together with equation CD determines the ionization and 
temperature of the gas for a given and N^. For this sxtuatron 
we use the total recombination coefficient in equation (1) , 
we take the Lya number density to be negligible. In this case only 
the ratio of these two densities is important. Some representative 
values are given in Table 1, in the column headed by "Extreme Dust" 
If the two-photon decay of the 2s state of hydrogen is the 
dominant cooling mechanism, then the rate at which energy is ra 
diated by the gss is given by 


(1 - fj^lNjjA'exp £-1.18 X lO^/^a] 


( 10 ) 


where A' is the transition probability for the two-photon decay, 
in this case the ionization and the temperature are determined by 
solving equations (1) and C7i simultaneously with the equation of 
expressions C8) and (10) . Results for this situation are also 

presented in Table !• 

once the ionization and temperature are determined, we can 
estimate the amount of ambient gas that is required for a parti- 
cular source: this estimate also depends on the cooling mechanism 
that dominates. Assuming that Lya is immediately degraded, the 
energy balance takes the form 


= Cl - fi)£iNH^7exc<''>^'’''o 
V 


( 11 ) 


Thus the required amount of gas is 

JE. 


NrV = 


- f >f.N <c (v) >hv 
^1 exc o 


( 12 ) 
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We can eaailjf seie tlii^ nuirihej: is fiox a, typical case. For 

the averag'e oEiaerved luBninosity in Broad. Baliner lines ^ we require J 
to be a few solar mass-es of suprathermals per year. So, if we 
choose E = 10® eV, f. = 0.1, and N„ = 10^®cm^^, about 100 solar 

O 1 ^ 

49 3 --5 

masses of qas is required in a volume of about 10 cm or 10 
cubic light years. We should note that if the optical reverberation 
process discussed in the next section is important, the value for 
J could be considerably reduced since the number of broad-component 
photons produced by each suprathermal proton could be greatly 
increased. 

If instead the two— photon decay is doing the cooling , the 
energy balance becomes 

= Cl - exp ^-1.18 X lo'^/T^hv^ (13) 


Again we can solve for the amount of gas needed: 


NjjV = 


Cl 


JE 


fi)A'hv^ 


exp 


11 . 


8 X 10 


(14) 


Using the same values as before and taking T^«?11,000°K, we find 

that a minimum of about 100 solar masses of ambient hydrogen is 

required in this case as well. 

Finally, it is possible to set a lower limit on the amount 

of gas needed irrespective of how the energy is ultimately radiated. 

This is because the solf— consistent calculation indicates that 

N /N < 10~®. If the relative density is larger than this, the 
s' II ^ 
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^hient becoj^e^ coiijplateljr ionized a,nd the ^toiuic processes 

are not efficient enough, to be of interest. If we coinbine this 
inequality" with equations C2I and C3I we obtain: 

/ E \ 

N„V > 5 X 10 ^^ — ^ I (15) 

^ f^Njj \^25kevy 

with the numerical values previously used/ this relation yields 
a lower limit of about lOM^ of ambient gas. 

If the physical conditions of the ambient gas are specified, 
we can begin to make estimates concerning dynamic stability . Here 
we examine the possibility that the force exerted on the gas by 
the suprathermal particles is balanced by the gravitational attrac- 
tion of a massive central body. We will consider the extreme case 
where all the suprathermals move in the same direction and have 
the same initial energy. 

Let us consider a cloud of ambient gas at a distance r from 

the source of the fast particles, which has an area A perpendicular 

to the motion of the suprathermals. Then the number of supra- 

thermals stopping in the cloud is — A. If the initial particle 

4irr^ 

energy is E^, and if all the particles stop in the cloud, then the 

J % 

rate of momentum transfer to the cloud is A(2m E ) . The 

4irr^ ^ 

most extreme case occurs when the cloud is just thick enough to 
stop the particles; this means that the volume of the cloud is 
A • R, where R is the stopping distance for the fast particles. 

For this case the rate of momentum transfer to unit volume of the 
ambient gas is 

J C^^p^o)^ 

4iTr^ 


R 


(16) 
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The gra,vtta,tional ;^Opce on unit volume of the gas due to the 
central source is' 


2 

r 


where M is the mass of the central object. Equating this force 
s 

to the momentum transfer rate gives us 


J 

4itR 




(17) 


So equilibrium obtains v;hen 



_ 1 _ 

Gmjj 


(18) 


The factor RN„ has previously been calculated for suprathermal 

protons in hydrogen gas (Ptak and Stoner 1973, figure 6), and it 

depends on and the ionized fraction of the gas f^. 

A few solar masses of protons per year are required to account 

for the observed Ha luminosity (the optical reverberation mechanism 

could reduce this) ? this makes J ^ lO^^sec If we assume that 

the initial energy of the protons is iMeV, and take the gas to be 

21—2 

10% ionized, then 10 cm . Inserting these values xn equa- 

tion (18) we obtain 


44 


gm 


5 X 10^® M . 

© 


< 10 
s 
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This is cjuite a largs it i.s o;^ the ordex o;^ magnitude of 

a spinar CHorrison and Cavaliere 1970i, W,e should note that this 
quantity depends more or less linearly on and that it goes 
like Furthermore, we have considered the most extreme 

case where the cloud of gas is just thick enough to stop all the 
particles. For the large ambient gas densities we have discussed, 
the value of R is quite small, and we should expect a typical gas 
cloud to be one or more orders of magnitude thicker than R. The 
value of Mg required by the equilibrixam condition would then be 
reduced by one or more orders of magnitude. 

IV. Radiation From a Suprathermal Atom 

A photon with frequency near the center of a Lyman or Balraer 
line will be repeatedly scattered by the ambient hydrogen atoms 
with small frequency changes due to the thermal motions of the 
atoms (which we assume are described by the temperature T^) . This 
scattering will both produce a distribution in frequency near the 
line centers (see Avery and House, 1968 and Auer 1968) and establish 
a quasithermal equilibrium between the distributions over the lower 
levels of the ambient atoms and the radiation field. As an appro- 
ximation, we assume that this equilibrium is established over the 
frequency range in which the gas is optically thick (t > 1) and 
that a photon scattered outside this frequency range always escapes. 
The frequency range in which the radiation field rises to the level 
of black body radiation then lies Av on either side of the line 
center , where 
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u 


Jin 




:a - f.)NjjLAX 




5 - 




(19) 


where U is the thermal doppler width in frequency = 


c 


^2kT ' 
e 


m 


H 


N„L is the hydrogen column density, A is the appropriate rate con- 
stant for absorption of the line, e is the energy of the lowest 
level for the transition relative to the ground state, and X is 
the wavelength of the photon. 

The suprathermal proton travelling through the ambient gas 
will spend a fraction of its time as a suprathermal atom, which 
can then interact by resonant absorption with the radiation field 
described above. In the rest frame of the suprathermal atom, the 
radiation field will not appear isotropic; only those photons tra- 
velling in directions approximately perpendicular to the supra- 
thermal atom's velocity will be available for resonant scattering. 
Thus, while an ambient hydrogen atom "sees" isotropic line radiation 
at the black body level, only a fraction AQ/4 it of the full solid 
angle of this radiation is available for resonant scattering with 
the suprathermal atom, and, at speed v, this is 


_ c Av 
4ir V 

o 


( 20 ) 


or unity when v is less than cAv/v^. At a speed of ac = 1 atomic 

40 

unit (2188 km/sec) , temperatures of about 10 K, an ambient density 
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of 10 cm , and for an ambient gas cloud 10 cm thick, is 

about 0.02 for the first few lines of the Lyman and Balmer series. 

At reasonably high temperatures and reasonably low densities, 
the fast atom will absorb and re-emit many line photons between 
atomic collisions; in this low density regime, the distribution 
over atomic states of the suprathermal atom will be determined by 
the size of AJ1/4 it and by the radiation temperature. At higher 
ambient gas densities certain inelastic collisions with the elec- 
trons, protons and atoms of the ambient gas begin to play a signi- 
ficant role. 

The rate at which the suprathermal atoms emit photons in a 
given line is the product of the probability the atom is in the 
appropriate upper state and an emission rate coefficient; thus, a 
calculation of the distribution over atomic states will yield pre- 
dictions for the relative intensities of the broad components of 
the various hydrogen emission lines. We have made such a calculation 
using a truncated set of eight hydrogen atomic states: Is, 2s, 

2p, 3s, 3p, 3d, all n = 4 states combined, and the ionized state. 

The relative probability the suprathermal hydrogen exists in each 
of these states is determined by equating the rates for populating 
and depopulating each one by the several radiative and collisional 
interactions at work. 

The radiative interactions include resonant absorption of Lya , 
Lv 3, Lyy, Ha and H3 from the radiation field as well as spontaneous 
emission of the same quanta- Collisional processes included are 
the excitations from the ground state, collisional mixing of the 
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2s and 2p states, population of all bound states by charge transfer 
from ambient atoms to suprathermal protons, and ionization from the 
ground state. Rates of other processes that could lead to transi- 
tions among these states, such as stimulated emission, were esti- 
mated and found negligible. Free parameters in the calculation 
were the radiation temperature, the ambient hydrogen number den- 
sity, the ambient ionization fraction, and the value of Af^/4TT at 
V = ac. The velocity-dependent cross sections used are generally 
those used for the previous computation of emission profiles (Ptak 
and Stoner 1973) and the collisional 2s-to-2p cross sections of 
Seaton (1955) . 

Table 2 gives some typical results of the calculation des- 
cribed above. The table presents quantities analogous to the "b's" 
commonly used to describe the distribution over states of hydrogen 
atoms in planetary nebulae. The quantity b(nl) is the ratio of 
the probability a suprathermal atom is in the state (nl) to the 
probability an ambient gas atom (in thermal equilibrium with a 
radiation field of temperature T^) is in the same state. If the 
suprathermal atom could interact resonantly with all the line 
photons rather than just those with momenta in the fractional 
solid angle AJ2/4 it, and if inelastic collisional processes were 
unimportant, then all of the b(nl) would be unity. 

Finally, Table 3 gives some typical results for the production 
of line photons by the suprathermal protons, obtained by integrating 
the production rate over the deceleration of the suprathermal proton 
from 6500 3cn/sec to about 500 km/sec. 
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An important result is that both the Lyman decrement and the 
ratio of the emission in Balmer lines to the emission in Lyman 
lines deoend strongly on the physical conditions in the ambient 
gas. This happens primarily because of the changing relative im- 
portance of atomic collisions and optical reverberation in this 
range of temperature and density. The importance of optical rever- 
beration is increased by either increasing the temperature (and 
thus the photon density in the radiation field) or decreasing the 
ambient gas density (decreasing the rate of atomic collisions) . 

T'Then optical reverberation dominates, the line emission rate 
per suprathermal is inversely proportional to density. At densi- 
ties of a few times lO^cm"^ and at T^ = lO'* K, a single suprathermal 
is capable of scattering 1 MeV of line radiation into the wings, 
so that this mechanism can easily remove all of the heat energy 
delivered to the gas by the suprathermals. This suggests the 
possibility that the temperature and density of the gas may adjust 
itself at and ^10 ^cm"^ to balance the heating and cooling 

effects of the suprathermals. 

In the limit of high density and low temperature, the emission 
line profiles will be the same as previously reported (Ptak and 
Stoner 1973) since the atomic collisions will the dominate the 
emission processes. As optical reverberation becomes important, 
the full width at zero intensity will remain the same since the 
characteristic speed at which the suprathermal proton spends a 
significant fraction of its time as a hydrogen atom will not change? 
the details of the emission profiles may change, however. 
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V. Discussion 

The picture that emerges is one in which a small volume of 
relatively high density gas is involved in the region where the 
broad emission lines originate. On the other hand, the ratio of 
the density of suprathermal particles to the density of gas part- 
icles is small; the upper limit is about lO”^. We find that the 
minimum amount of ambient gas which is required is moderate, about 
100 for an ambient gas density of 10^*^cm~^. 

Under these conditions, we can construct a self-consistent 
picture in which the temperature and ionization of the ambient gas 
are determined by the flux of suprathermal protons, and where the 
hydrogen emission excited in the gas is trapped and degraded before 
r it can escape. By finding reasonable physical conditions in a 
nuclear gas, consistent with the presence of suprathermal particles 
and with the removal of the "narrow component" of Balmer radiation, 
we have removed two major objections to the suprathermal ion hypo- 
thesis for the origin of the broad line emission. The physical 
conditions of the gas will depend somewhat on the mechanism which 
converts the kinetic energy of the suprathermals into, radiation. 
Detailed self-consistent calculations have been performed by 
Kimmer (1974) in two limiting cases: 1) the Lya photons excited 

in the gas are immediately destroyed, 2) the energy escapes from 
the gas only by way of the two-photon decay of the hydrogen 2s state. 
His results are qualitatively the same as those presented in Table 
1, for which we have assumed that the gas and the suprathermal 
protons are everywhere the same. 
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The first case provides us with a lower limit for the temper- 
ature and ionization of the gas^ and the second case gives us an 
upper limit for the electron temperature and ionization fraction. 
Intermediate cases could be provided by varying amounts of dust 
in the gas absorbing Lya photons and reradiating the energy in the 

infrared . 

Another possible cooling mechanism which lies between the two 
extreme cases is the optical reverberation from suprathermal hydro- 
gen atoms. We have not yet calculated the physical conditions of 
the gas in a self-consistent way with this process dominating the 
cooling of the gas. This does not seem crucial since the two 
limiting cases do not produce radically different results. What 
we have calculated are the population of states for the suprathermal 
atoms as a function of the temperature and density of the ambient 
gas . 

With these population ratios available, it is a simple matter 
to calculate the intensity ratios for the broad components of the 
hydrogen emission lines. The calculations we have done so far 
determine only the Lyman decrement reliably (since the angular 
momentum sublevels of n = 4 are not treated separately) , and these 
results appear in Tables 2 and 3 . This intensity ratio has not 
yet been measured for any Seyfert galaxy, of course. However, we 
can see from the table that the calculated ratio is very sensitive 
to the density and temperature of the ambient gas, and we can expect 
that the Balmer decrement calculated in this way will also be a 
sensitive indicator of the physical conditions in the gas. We are 
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presently working to extend the calculation in this direction. 

The dynamic picture is still rather unclear. We have seen 
that a kind of force balance is possible if a very massive central 
object binds the ambient gas gravitationally. How massive this 
object must be depends on the geometry and on the sorts of res- 
trictions one wishes to place on the size of the emitting region. 
rJeometrical factors have been suppressed in the model presented 
in this paper. However, it is possible to make some statements 
of this type based on our calculations. 

We found that for an ambient gas density of about lO^^cm ^ 
and an ionized fraction of 0.1, the required volume of gas is 
about lO^^cm^. So if the emitting region is to be confined to 
within a sphere of radius one light year, the minimum value of 
the filling factor is about 2 x lO”^. This small value implies 
that even with a small size for the emission region, a good deal 
more ambient gas may be present than is needed, or, in other words, 
the qas clouds may be considerably thicker than the stopping dis- 
tance of the suprathermal particles. The necessary mass for a 
central object to produce some quasi-static equilibrium is then 
correspondingly reduced. 

Another way of viewing the problem is to ignore gravity and 
to consider only inertia. With the conditions described above, 
the minimum value of ambient gaS in the emission region is about 
100 M^. If this is all the gas present, and if about 1 of 
fast protons with parallel initial velocities stop in the gas each 
year, then the system is dynamically unstable with a time scale on 
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the order of 100 years. However, if there is much more gas present 
than the minimum required, then the time scale is correspondingly 
increased . 

An important feature of the self-consistent calculation is 
that when the relative density is not so large as to make the gas 
completely ionized, then for a wide range of conditions, the 
ambient gas is quite neutral; the ionized fraction is like 0.1 
or less. This means that when the atomic collision mechanism is 
able to operate, it does so at close to the best efficiency pos- 
sible. It also means that a moderate amount of gas can have a 
large optical thickness in Lyman and Balmer lines, and thus pro- 
vide the necessary trapping of the line photons excited in the 
ambient gas. 

We are in the process of calculating detailed line profiles 
and the Balmer decrement for the case where optical reverberation 
dominates the emission from the suprathermal particles. The re- 
sults of this calculation will be presented in a subsequent paper. 
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Table 2 


The distribution over the lowest few excited states of a suprathermal hydrogen atom with 
velocity v = ac in an optically thick hydrogen gas. In all cases the gas is 5 percent ionized 
and the quantity AJ^/4ir = 0.02. See text for the definition of the b(nl) . 
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Table 3 


The total line radiation from a suprathermal hydrogen stopping in an optically thick hydrogen 
gas under different physical conditions. In all cases the quantity in equation (20) is taken 
to be 0.02 for v = ac. 
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